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Abstract 

We have studied the dynamics of spreading of viscous non-volatile fluids on 
surfaces by Monte Carlo simulations of solid-on-solid (SOS) models. We 
have concentrated on the complete wetting regime, with surface diffusion 
barriers neglected for simplicity. First, we have performed simulations for the 
standard SOS model. Formation of a single precursor layer, and a density 
profile with a spherical cap shaped center surrounded by Gaussian tails can 
be reproduced with this model. Dynamical layering, however, only occurs 
with a very strongly attractive van der Waals type of substrate potential. 
To more realistically describe the spreading of viscous liquid droplets, we 
introduce a modified SOS model. In the new model, tendency for dynamical 
layering and the effect of the surface potential are in part embedded into 
the dynamics of the model. This allows a relatively simple description of 
the spreading under different conditions, with a temperature like parameter 
which strongly influences the droplet morphologies. Both rounded droplet 
shapes and dynamical layering can easily be reproduced with the model. 
Furthermore, the precursor width increases proportional to the square root 
of time, in accordance with experimental observations. 
PACS numbers: 68.10.Gw, 05.70.Ln, 61.20.Ja. 
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1 Introduction 



Surface processes such as wetting coating and adhesion depend strongly 
on phenomena occuring in the first few layers on top of the substrate. Recent 
progress in optical techniques [2-8] has focussed attention to the dynamics 
of spreading of tiny, viscous nonvolatile droplets on solid surfaces. With 
ellipsometry, measurements can now be done on molecular scales. They have 
revealed fascinating phenomena previously unaccounted for in macroscopic 
droplets, such as dynamical layering and diffusive spreading of the precursor 
layer EE^E- 

In the experiments, the most common liquids have been tetrakis (2-ethylhex- 
oxy)-silane, squalane, and polydimethylsiloxane (PDMS). It has been shown 
1^ that the thickness profiles of tetrakis and PDMS droplets on a silicon wafer 
exhibit strikingly different shapes under spreading. On a "low energy" surface 
(with a relatively low solid-vacuum surface tension fl]]) tetrakis produces 
clear dynamical layering (Fig. 1(a)), i.e. stepped shapes of exactly one 
molecular layer in thickness. In contrast, on the same substrate the spreading 
of PDMS proceeds by a fast evolving precursor of about one molecular layer 
in thickness (Fig. 1(b)). 

On the other hand, the same liquid on different surfaces has been shown to 
form distinct morphologies |Q . On a "high energy" surface (with a relatively 
high surface tension and adsorption energy |jl|) PDMS developes dynamical 
layering, too, in analogy to tetrakis. On a low energy surface the formation of 
a single precursor layer with a central cap is observed . Other experiments 
show P, |] that PDMS droplets form stepped, dynamically layered profiles 
on two different silicon wafers (Fig. 1(c)). In addition, a recent experiment 
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of PDMS spreading on a silver substrate shows the formation of a profile 
with spherical cap shaped center and Gaussian tails at late times (Fig. 1(d)). 



On the theoretical side progress has been more moderate. The continuum 



hydrodynamic approach 0, [TO, O, while being able to qualitatively ex- 



plain the precursor film and its eventual diffusive (~ t^/^) growth |Tl|] is 
simply not valid for the measurements on scale of Angstroms. The first the- 
oretical model to successfully account for dynamical layering was given by 
Abraham et al. [|T^. They used the two dimensional horizontal solid-on- 
solid (SOS) model, which describes coarse-grained layers of fluid emanating 
from a droplet "reservoir". Dynamical layering was observed to occur with 
an attractive van der Waals type of surface potential in the complete wetting 
regime [0, |13|. The model also predicts diffusive behaviour for the bulk 
layers of the droplet on a high energy surface, but linear time dependence 
for the precursor layer. 

A different approach was recently taken by de Gennes and Cazabat [Q, 
who proposed an analytic model of an incompressible, stratified droplet. 
The stepped morphology of the spreading droplet is postulated a priori, and 
spreading is assumed to proceed by a longitudinal flow of the molecular lay- 
ers, and a transverse permeation flow between them. The model introduces 
an effective Hamaker constant, which has been fitted to experimental data 
on stepped PDMS droplets. 

On a more microscopic level, molecular dynamics simulations have been per- 
formed WB, W^, O. Yang et al. |]TB[ studied droplets of Lennard- Jones (LJ) 



particles and dimers on corrugated surfaces, surrounded by a vapor phase. 
They observed terraced droplet shapes in the complete wetting regime, but 



obtained logarithmically slow growth for the layers. In contrast, Nieminen 
et al. [|1^ were able to qualitatively observe a crossover from "adiabatic" to 
"diffusive" spreading for LJ droplets on a smooth surface. To account for 
the chainlike molecular structure of some of the liquids, Nieminen and Ala- 
Nissila |jl8| recently carried out a systematic study of droplets of LJ solvent 



particles and flexible oligomers on smooth surfaces. By varying the strength 
of the surface potential, they were able to qualitatively conclude that the 
actual molecular structure of the liquid can play an important role in the 



formation of the density profiles, even beneath the entanglement regime [|T^. 
Namely, relatively weaker substrate interactions and longer chains favour 
more rounded droplet shapes, due to the mixing of the molecular layers un- 
der spreading. On the other hand, a very strong surface attraction promotes 
layer separation and leads towards stepped shapes. In the latter case, the 



spreading process qualitatively proceeds in analogy to the ideas of Ref. WA 



with layers of fluid "leaking" at the edges of the steps. They also observed a 
crossover from "adiabatic" to "diffusive" behaviour of the precursor film and 
developed a scaling form for its time dependence. These results are in good 
qualitative agreement with experiments; however, droplets consisting of just 
a few thousand molecules were simulated, while experimental droplets are 
microscopic in the vertical direction only. 

In the present work we have taken to approach the problem of droplet spread- 
ing on a more coarse-grained level, in order to further clarify the role of the 
effective viscosity of the liquid, and the surface potential in the droplet mor- 
phologies. We have carried out systematic Monte Carlo simulations with 
two different SOS models to study the dynamics of the spreading process. In 
these models, the droplet is assumed to consist of effective, coarse-grained 
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blocks of molecules |T^. The size of these blocks can be highly anisotropic 
between the vertical and horizontal directions, depending on the details of 
the interactions. All the simulations have been done in the complete wetting 
regime in the canonical ensemble, where the final state of the droplet is a 
surface monolayer. No evaporation from the surface is allowed, and in this 
work the surface diffusion barriers have been neglected for simplicity. 



The first model that we have studied is the standard vertical SOS model [20 



21| . Besides computing density profiles of the droplets, we have also followed 
the time evolution of the width of the precursor film r{t), and the droplet 
height. Our results demonstrate that with a strongly attractive short range 
substrate potential, only a single precursor layer evolves, with the center of 
the droplet assuming a spherical cap shape rather than a Gaussian profile. 
At later times, a surface monolayer forms. These results are analogous to 
the observations in Ref. |2T[]. In addition, formation of a spherical cap with 
Gaussian tails occurs with relatively weak attractive short range substrate 
potentials. To obtain stratified droplet shapes corresponding to dynamical 
layering, however, an extremely strong long range van der Waals type of 
surface potential must be used. In addition, the behaviour of the precursor 
film was not observed to follow the correct dynamical behaviour. This is due 
to the somewhat unrealistic dynamics of the model. 

To better take into account the fluid-hke nature of the spreading droplets, we 
have developed a modified SOS model in which the tendency for dynamical 
layering and the effect of the surface potential have been partially embedded 
into the dynamics of the model, in analogy with some of the ideas presented 
by de Gennes and Cazabat In its simplest form, the model enables a de- 



scription of the spreading process with a single temperature like parameter. 



When varied, this parameter strongly influences the droplet morphologies 
and reproduces both rounded and stepped droplet shapes. In particular, dy- 
namical layering can be reproduced without explicitly including a long range 
van der Waals substrate attraction. When such an attraction is added the 
range of layering is further enhanced, as expected. Our quantitative results 
also show that the precursor width increases proportional to the square root 
of time. A preliminary account of these results has been published in Ref. 
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2 Standard Solid— on— Solid Model of Spread- 
ing Dynamics 

2.1 Definition of the Model 

The standard solid-on-solid (SOS) model used in our simulations is deflned 
by the following Hamiltonian [^: 

i,j i=2 i 

where in the flrst term r < is the nearest neighbour interaction parameter 
and the summation goes over the four nearest neighbours of the column i 
whose height is hi. The heights, which measure the average concentration of 



particles and obey a vertical SOS restriction are integers. The second 
term on the right hand side of Eq. (1) describes the vertical interaction 
between the effective molecules. V{hi) denotes the strength of an attractive, 
height dependent substrate potential at site i. V{hi) is deflned as 



V{hi) = ^ + Bh,5h^,u (2) 
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where A < is an effective Hamaker constant, and B < describes the 
interaction between the substrate and the fluid. 

In the model, the effective particles follow standard Monte Carlo (MC) 
Metropolis dynamics with transition probability 

P{hi ^ h'i) = min{l, e-^^/'^^^}, (3) 

where AH is the energy difference between the final {h'^) and initial {hi) 
states (with the SOS restriction obeyed) . Time is measured in units of MC 
steps per site (MCS/s). The initial configuration of the droplet is usually 
chosen to be a three-dimensional cube, which after a transient time assumes 
its characteristic shape. During the spreading stage, we calculate the density 
profile of the droplet (in analogy to the experiments), the time dependence 
of the precursor width r{t) and the height of the droplet. All the simulations 
have been done without allowing evaporation from the surface; thus the final 
state is a surface monolayer. In this work the activation energy for surface 
diffusion of isolated particles is neglected, and thus the model may not be 
applicable to some high energy surfaces. 

2.2 Results for the Standard Model 

First, we study the spreading dynamics without the van der Waals potential 
by setting A — in Eq. (2). In this case, the results are not very sensitive 
to the droplet size, and in most cases allxllxlO droplet was used, where 
ten is the initial height. We also fix r —0.05 in this section. With a 
relatively strong surface attraction, only a single precursor film is formed, 
with a rounded central cap. At the edge of the film the migration of the 
effective molecules causes the formation of a monolayer of separate particles. 
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This happens typically in the range \B /t\ = 2.0—4.0, and \kBT /t\ = 0.6—1.0. 
These results are in complete agreement with the qualitative observations of 
Ref. Additionally, we find that the time development of the precursor 

film approximately follows r(t) ~ t"'^^, in contrast to t^^'^. If we then lower 
the surface attraction and set B = t, and \kBT/T\ = 0.6, the precursor film 
vanishes and the droplets assume a rounded, spherical shape with Gaussian 
tails as shown in Fig. 2. 

Next, we study the infiuence of the long range van der Waals type potential. 
For a moderately large A, the results are virtually indistinguishable from 
those above. Namely, for A = B = —0.15, a single precursor film forms while 
for A = B = T the film disappears {\kBT/T\ = 0.6). However, for a very large 
surface attraction A = B = —1.0, the droplets show indications of separation 
of layers. At IksT/rl = 1.0 shown in Fig. 3(a), a small central cap is formed 
at submonolayer heights. When temperature is further lowered to 0.6, three 
separate layers can be observed (Fig. 3(b)). We have calculated r{t) ~ for 
the lower temperature case, and find a ~ 0.17 for both A = B = —0.15 and 
— 1.0. We also find that at very late times, the height of the submonolayer 
film decreases approximately as t~^. This describes the decrease of the almost 
two dimensional island, which leaks into a surface monolayer. This monolayer 
which spreads diffusively consists of all the particles, which are two or more 
lattice sites apart from the cluster. Thus the precursor film width decreases 
and the areal density of the monolayer correspondingly increases. We note 
that due to the lack of surface diffusion activation energy, the case of larger 
surface attraction enhances the rate of spreading, which is opposite to what 
happens on high energy surfaces |TB| . 



From the results presented above, it is clear that the standard SOS model 



can produce dynamical layering only for exteremely large values of the long 
range surface attraction in the model. In fact, although standard SOS mod- 
els have been very useful in describing solid surfaces and their roughening 
20[| , and even growth of surface layers it is clear that they are less re- 



alistic for fluid spreading. For example, the diffusion length of the effective 
particles is not well-defined, due to the fact that very large height differences 
between nearest neighbour are allowed. Second, nonmonotonic variations in 
the density profiles are possible, although not significant for the parameter 
range used in the present work. Due to these limitations, in the next sec- 
tion we shall introduce a modified SOS model which better accounts for the 
fluid-like nature of the spreading droplets. 

3 Modified Solid— on— Solid Model of Spread- 
ing Dynamics 

3.1 Definition of the Model 

The basic idea behind the modified SOS model is to more realistically de- 
scribe the viscous fiow of the droplets. We have done this by modifying the 
dynamics of the standard SOS model defined by Eqs. (l)-(2) in the following 
way. First, the effective molecules are not allowed to climb over height barri- 
ers, or jump down over terrace edges. For example, two individual particles 
located initially at different sites on the surface can never form a column of 
height hi = 2. Second, when the difference between the neighbouring initial 
and final sites is {hi — h'^ > 2, a transient non-SOS excitation (a hole or 
vacancy) can be created into the initial column by a particle which jumps to 
the neighbouring position with a probability 
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P = mm{Q{hi - hj - 1/2), e-^^*/'^^^} (4), 

where Q{x) is the Heavyside step function. The energy difference AHf is 
calculated for this transient state. The hole left behind is immediately filled 
by the column of particles above it, which are all lowered by a unit step. 
This means that AHt > AH (as calculated for the final state). We note that 
the effective dynamics used here does not satisfy microscopic reversibility at 
all times. 

The motivation for the unusual dynamics comes from the physics of dissipa- 
tive spreading of viscous fluids on attractive surfaces. Microscopic calcula- 



tions have revealed |T8[ that at least in the layered state, the droplets spread 
by the viscous flow of rather well separated layers, which "leak" on the surface 
at the terrace edges. Furthermore, the analytic model for stratified droplets 



by de Gennes and Cazabat |TJ] assumes, that the horizontal layers of fluid 
flow with friction, and that there is a narrow region of permeation flow at 
the terrace edges. In our model, this is roughly described by the creation 
of the transient hole at the step edges. Since this is temperature controlled 
(cf. Eq. (4)), we expect that the model be able to describe droplet spreading 
under a variety of different conditions. 

3.2 Results 

First, we present results for A = B = r = 0, which corresponds to an infinite 
temperature. The effective dynamics in the model renders even this sim- 
plest case nontrivial. As seen in Fig. 4, at early and intermediate stages of 
spreading the droplet assumes a non-Gaussian, spherical shape with Gaus- 
sian tails surrounding the center of the droplet. At early times the spreading 
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is relatively fast, but slows down at later times as the droplet density pro- 
file approaches a Gaussian shape at submonolayer heights. In contrast, we 
note that an infinite temperature for the standard SOS model always leads 
to Gaussian shaped profiles. At late times, the width of the precursor film 
approaches the expected t^^'^ behaviour. 

Next, we set A — B — T — —0.05 and let the temperature ratio ksT/T 
vary. For the lowest temperature of 0.6, dynamical layering of the density 
profiles is immediately evident in Fig. 5(a). For larger droplets, up to five 
or six layers can easily be detected (Fig. 5(b)). In this case the droplet 
maintains a rather compact shape and the surface monolayer formation is 
relatively slow. At later stages (not shown in Fig. 5) a flat proflle is formed 
with a tiny central cap at the center. This cap disappears at approximately 
after 6000 MCS/s. The precursor width follows ^0-50(3) ^Iso studied the 
spreading of an initially ridge-shapcd droplet of size 120 x 11 x 10, where the 
precursor width follows t^^'^^ rather than t^/'^. In this case the spreading is 
highly anisotropic and rather different from the three-dimensional droplet. 

When temperature is raised, the layer formation weakens and droplets be- 
come more spherical. For l/csT/rl = 0.8 (Fig. 6) and 1.0, the precursor film 
grows as t^-^^^'^) and ^^■^'''(s)^ respectively. The surface monolayer formation 
due to enhanced thermal fiuctuations is faster, and for the higher tempera- 
ture the whole precursor film tends to break down. In Fig. 7 we show the 
time development of the droplet height at the three different temperatures. 
The fast spreading rate, and consequently the enhanced edge flow leading 
to more rounded shapes at higher temperatures is clearly reflected in the 
results. 

Next we study the spreading with a long range van der Waals type of surface 
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potential with 1^4/5 1 = 1. In this case simulations have been performed for 
B = -0.025, B = -0.05 = r, B = -0.07, and B = -0.2. The tempera- 
ture was set to IksT/rl = 0.6. In our model with modified dynamics, the 
inclusion of the van der Waals energy enhances the tendency for dynamical 
layering, as demonstrated in Fig. 8. Furthermore, the results indicate clearly 
that with B = —0.2 the precursor films breaks down, giving for the time de- 
velopment of the precursor film width the behaviour ~ With the other 
choices, however, the time dependence of the precursor width follows t^^^ 
rather accurately. For B = —0.2 spreading becomes unrealistically fast due 
to the lack of the surface diffusion activation barriers, which have been set 
to zero. As mentioned earlier, we do not expect the model to be applicable 
to high energy surfaces if the diffusion barriers control the dynamics. 

4 Discussion of the Results and Comparison 
with Experiments 

Next, we analyze the density profiles and make comparisons with relevant 
experiments [2-8], whenever possible. Due to the coarse-grained nature of 
the SOS models quantitative comparisons may not be meaningful, however. 
The aim here in part has been to eludicate the physical processes behind 
different droplet morphologies, and thus we will also discuss our results in 



connection to other coarse-grained models |12, 13, 14, 21 1, and microscopic 



simulations iTGl, ITTl 



First, we discuss results for the standard SOS model. The model quali- 
tatively reproduces single precursor film formation with a relatively strong 
attractive short range substrate potential, while with a weaker potential, 
rounded droplet shapes with Gaussian tails appear. This is in accord with 
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the notion that increased surface attraction should enhance the precursor 
film ^, |12|, [1^. The precursor film was also qualitatively observed in 



the simulations of Ref . . Also, with the inclusion of a long range van der 
Waals potential, the droplets exhibited tendency towards dynamical layering. 
This is in accordance with the horizontal SOS model results of Refs. iO, ITHl, 



and qualitatively with the microscopic picture as well |18[. Although 
the standard SOS model can reproduce density profiles resembling experi- 
mentally observed ones (Fig. 3(a) and Fig. 3(b) @, P)) the strength of 
the long range potential must be considered unrealistically large as compared 
with the particle-particle interactions on the surface level. Most importantly, 
the time dependence for the width of the precursor film does not follow the 
expected t^^^ behavior, which reflects the unrealistic nature of the dynamics 
in the model. 

In contrast, the modified SOS model yields results which seem much more 
realistic. At low temperatures, stepped droplet shapes are obtained (Figs. 
5(a) and (b)). On the other hand, when particle motion from one layer to an- 
other is more likely, rounded droplet shapes occur (Figs. 4 and 6). This is in 
qualitative accord with conclusions from coarse-grained models [l^, |I4 



and microscopic calculations ||17| , |T8[. Furthermore, the inclusion of a van 
der Waals potential further enhances the range of layering, as expected. The 
results are also in qualitative agreement with experiments 0, ^ ^ where 
increasing the viscosity of the liquid enhances layering, and layering is also 
more prominent on high energy surfaces. Most importantly, the model re- 
produces the expected t^/"^ dependence for the width of the precursor layer. 
However, it should be noted again that since the model neglects diffusion 
barriers, our results may not be applicable for all cases. 
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The remarkable feature of the modified SOS model is that by just varying 
a single temperature like parameter, droplet morphologies from rounded to 
stepped shapes can be easily reproduced. This can be understood based 
on the physical analogy behind the modified dynamics [Q. High effective 



temperatures in the model correspond to enhanced flow at the edges, and 
reduced viscocity and surface interactions allowing transient hole formation 
with relative ease. On the other hand, lowering the temperature suppresses 
this hole formation which is analogous to suppressing interlayer flow in the 
phenomenological model of de Gennes and Cazabat In fact, our model 



gives further support to the physical ideas behind their model regarding 
dynamical layering. 

5 Summary and Conclusions 

To summarize, in this work we have tried to unravel physical processes be- 
hind different morphologies observed in the spreading of tiny liquid droplets 
on surfaces. The present approach of using coarse-grained descriptions of 
the fluid is complementary to other approaches, and can indeed give a lot of 
physical insight to the problem. In particular, by introducing a new, more 
realistic coarse-grained lattice SOS model of fluid spreading we have demon- 
strated how the interplay between the internal viscosity of the fluid and the 
nature of the surface attraction can produce a number of different droplet 
shapes, in qualitative agreement with experiments. Our results also give fur- 
ther support to ideas behind the analytic models, in particular that of Ref. 



n| regarding the spreading of terraced, stratified droplets. They are also 
in agreement with microscopic studies of Refs. [0, 0. However, the ex- 
perimental situation remains complicated and should not be oversimplified 



13 



- there exist very few attempts to systematically study the effects of surface 
attraction, microscopic structure of the liquid, finite size effects and other 
factors, which can play an important role in spreading dynamics. We hope 
that the present work inspires more systematic studies in these directions. 
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Figure captions 



Fig. 1. (a) Experimental thickness profiles of tetrakis on a "low energy" 
silicon surface 0. Dynamical layering occurs at late times, (b) Thickness 
profiles of PDMS on the same surface (the "Mexican hat" shape) 0. For even 
smaller droplets, the central cap becomes more distinct and a Gaussian shape 
is assumed at latest times 0. (c) Thickness profiles of PDMS on a silicon 
wafer covered by a grafted layer of trimethyl groups. This surface behaves 
as a "low energy" surface with strong dynamical layering occuring. (d) 
Thickness profiles of PDMS on a silver surface after 50, 80, 150, and 360 
minutes The submonolayer profiles can be fitted by a spherical cap, with 
Gaussian convolution at the edges. 

Fig. 2. Density profile for the standard SOS model, with A = 0, B = t, and 
IksT/rl = 0.6. The droplet assumes a rounded profile, with no precursor 
film (cf. Fig. 1(d)). The initial height of the droplet is 41 units, and the 
result is an average over 50 runs after 25 000 MCS/s. The solid line indicates 
a Gaussian fit. 

Fig. 3. Spreading droplet profiles for the standard SOS model with A = B = 
—1.0, and (a) I^bT/tI = 1. At submonolayer heights, a central cap forms. 
The curves correspond to 800, 1200, 1400 and 1600 MCS/s, and are averages 
over 500 runs. The initial size of the droplet is 11 x 11 x 10. (b) Results for 
l/c^T/rl = 0.6. Dynamical layering tends to occur, with three layers visible. 
The curves correspond to 1600, 2200, and 3800 MCS/s, with averages over 
300-400 runs. 

Fig. 4. The density profile of a spreading droplet for the modified SOS model 
at an infinite temperature, after 650 MCS/s. A distinctly non-Gaussian, 
spherical shape occurs above monolayer heights, as evidenced by a Gaussian 
fit shown with a solid line. The initial height of the droplet is 41 units, and 
the profile is an average over 500 runs. 
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Fig. 5. Density profiles for the modified SOS model at |A;bT/t| = 0.6. 
Dynamical layering of the profiles occurs after an initial transient, (a) Results 
for a 11 X 11 X 10 droplet, with curves corresponding to 1800, 2800, and 3800 
MCS/s. (b) Results for a 11 x 11 x 40 droplet, with curves at 5800, 6800, 
and 8000 MCS/s. All profiles are averages over 500 runs. 

Fig. 6. Density profiles for the modified SOS model at \kBT/r\ = 0.8. A 
shoulder and a central cap are clearly seen, with a spherical cap shaped center 
and a Gaussian foot at late times. The curves correspond to 2000, 2200, and 
2600 MCS/s. Results are averages over 500 runs. 

Fig. 7. The height of the droplets for the modified SOS model at the three 
different temperatures studied, as a function of time. Filled squares corre- 
spond to I/cbT/tI = 0.6, open circles to 0.8, and filled circles to 1.0. Breaks in 
the curves correspond to separate layers. While the spreading rate strongly 
increases as a function of temperature, the dynamical layering considerably 
weakens and droplets become more rounded. 

Fig. 8. Density profiles for the modified SOS model, with |A/-B| = 1, 
IksT/rl = 0.6, and B = —0.2. A stronger van der Waals potential clearly 
enhances dynamical layering (cf. Fig. 5). See text for details. 
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